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Abstract: We compute the next-to-next-to-leading order (NNLO) QCD corrections to 
the first five moments of six event shape variables related to three-particle final states in 
electron-positron annihilation; the thrust, the heavy jet mass, the C-parameter, the wide 
and total jet broadenings and the three-to- two-jet transition parameter in the Durham 
algorithm Y%. The NNLO corrections to the first moment are moderate for all event shapes, 
while the renormalisation scale dependence of the theoretical prediction is substantially 
reduced compared to the previously existing NLO results. From a comparison with data 
from JADE and OPAL, we observe that the energy dependence of the moments of the 
wide jet broadening and I3 can be largely explained without any non-perturbative power 
corrections, while the other observables exhibit a clear need for power-like contributions at 
low centre-of-mass energy. 
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1. Introduction 



Event shape variables in e + e~ annihilation provide an ideal testing ground to study Quan- 
tum Chromodynamics (QCD) and have been measured and studied extensively in the last 
two decades. In particular, event shape variables are interesting for studying the interplay 
between perturbative and non-perturbative dynamics. Apart from distributions of these 
observables, to which the NNLO corrections recently have become available [1—6], one can 
also study mean values and higher moments. The nth moment of an event shape observable 
y is defined by 



where y max is the kinematically allowed upper limit of the observable. Moments were 
measured for a variety of different event shape variables in the past. The most common 
observables y of three-jet type are: thrust T [7] (where moments of y = (1 — T) are taken), 
the heavy jet mass p = Mjj/s [8], the C-parameter [9], the wide and total jet broadenings 
By/ and Bt [10], and the three-to-two-jet transition parameter in the Durham algorithm 
I3 [11]. Definitions for all observables are given in, for example, Ref. [1]. Moments with 
n > 1 have been measured by several experiments, most extensively by JADE [12] and 
OPAL [13], but also by DELPHI [14] and L3 [15]. A combined analysis of JADE and OPAL 
results has been performed in Ref. [16]. 

As the calculation of moments involves an integration over the full phase space, they 
offer a way of comparing to data which is complementary to the use of distributions, where 
in general cuts on certain kinematic regions are applied. Furthermore, the two extreme 
kinematic limits - two-jet-like events and multi-jet-like events - enter with different weights 
in each moment: the higher the order n of the moment, the more it becomes sensitive to 
the multi-jet region. Therefore it is particularly interesting to study the NNLO corrections 
to higher moments of event shapes, as these corrections should offer a better description 
of the multi-jet region due to the inclusion of additional radiation at parton level. 

Moments of event shape variables can be divided into a perturbative and a non- 
perturbative part, 



where the non-perturbative part accounts for hadronisation effects. The non-perturbative 
part is suppressed by powers of X p /Q p (p > 1), where Q = y/s is the centre of mass 
energy and Ai is of the order of Aqcd- The functional form of X p has been discussed quite 
extensively in the literature, but as this parameter is closely linked to non-perturbative 
effects, it cannot be fully derived from first principles. 

The power corrections can be related to infrared renormalons in the perturbative QCD 
expansion for the event shape variable [17-24]. The analysis of infrared renormalon ambi- 
guities suggests power corrections of the form \ p /Q p , but cannot make unique predictions 
for \ p : it is only the sum of perturbative and non-perturbative contributions in (1.2) that 
becomes well-defined [25]. Different ways to regularize the IR renormalon singularities have 
been worked out in the literature. One approach is to introduce an IR cutoff u/ and to 




(1.1) 



(y n ) = (y n ) P t + (y n ) 



(1.2) 



- 1 - 



replace the strong coupling constant below the scale pi by an effective coupling such that 
the integral of the coupling below pj has a finite value [23, 24] 

1 f^ 1 

— / dQa cS (Q 2 ) = a (pi) . (1.3) 
Vi Jo 

This dispersive model for the strong coupling leads to a shift in the distributions 

<"> 

where the numerical factor a y depends on the event shape, while V is believed to be 
universal and scales with the CMS energy like pi/Q. We refer the interested reader to 
Refs. [24, 26, 27] for further details. 

Another ansatz is that of Korchemsky et al. [28,29] who suggest a shape function valid 
equally for thrust, C-parameter and heavy jet mass p, which is independent of the energy 
scale and takes into account the energy flow into two hemispheres of the final state. In the 
region y » Aqcd/Q, their predictions for the mean values of 1 — T and p coincide with 
those of renormalon based models, i.e. the leading power corrections are parametrised by 
a single non-perturbative scale Ai. In the case of the C-parameter however, they find an 
additional contribution modifying the size of the 1/Q correction. This formalism also allows 
one to study the region y ~ Aqcd/Q, where power corrections of the form (Aqcd/(Qu)) p 
for arbitrary p become important. 

Gardi et al. [30-33] make an analysis using the formalism of the so-called "dressed 
gluon exponentiation" [32], which resums the renormalon contribution to the Sudakov 
exponent, taking into account factorially enhanced subleading logarithms. Within this 
formalism, power law corrections beyond the leading term are predicted which probe the 
region y ~ Aqcd/Q- In particular, this formalism predicts that power corrections of the 
form X2/Q 2 are suppressed, which is also indicated by the data. Ref. [33] contains a 
detailed analysis of the power corrections for the thrust and heavy jet mass distributions, 
taking also hadron mass effects into account and deducing from renormalon analysis that 
the non-perturbative corrections to both distributions can be described by a single shape 
function. In Ref. [31], characteristic functions describing the leading and subleading power 
corrections to the first four moments of the thrust variable are derived. 

In Ref. [34], a renormalisation group improved (RGI) treatment of the first moment 
of event shape observables is suggested and compared to data. A simultaneous fit of three 
parameters - Aqcd, a parameter p2 describing uncalculated perturbative higher order 
corrections and a parameter Kq controlling the power corrections - showed for the cases of 
(1 — T) and (p) that the data are consistent with Kq = within the RGI approach without 
leading to unreasonably large values of p2, but do not allow definite conclusions. 

Overviews on the various approaches to estimate non-perturbative corrections can be 
found in Refs. [25,35]. Recently, significant progress in the theoretical description of power 
corrections has also been made using Soft-Collinear Effetive Theory (SCET [36-38]), see 
e.g. Ref. [39]. 

An extensive comparison of e + e~ data on events shapes with NLO+NLLA QCD pre- 
dictions, supplemented by power corrections within the approach of [23, 24], has been 
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performed in [40] . A determination of a s and uq from fits to distributions and mean values 
of five event shapes has been carried out, and indications were found that uncalculated 
higher orders may contribute significantly to the non-perturbative corrections. A determi- 
nation of a s based on moments of event shapes calculated up to NLO has recently been 
reported in [41]. 

Based on the NNLO corrections published in [1], several new determinations of the 
strong coupling constant have become available recently. In Ref. [42], fits of six event 
shape distributions to ALEPH data have been made. In [43], the same perturbative results 
have been fitted to JADE data, including also matching [44] onto the resummation of large 
logarithms [45] in the next-to- leading log approximation (NLLA). 

A very recent study of non-perturbative contributions to the thrust distribution, based 
on the NNLO corrections published in [1] and including NLLA resummation, as well as 
power corrections within the low-scale effective coupling formalism [23,24], can be found 
in [46]. This work also provides a determination of a s and «o based on fits to a wealth of 
experimental data for the thrust distribution and confirms the conjecture of [40] that the 
explicit 1/Q corrections become smaller as higher perturbative orders are included. This 
corroborates the hypothesis that the divergent renormalon contribution can be regularized 
by modifying the strong coupling at low scales, leading to the effective coupling used in 
eq. (1.3). 

In this paper, we provide the NNLO corrections to the first five moments of the event 
shapes thrust T, the normalised heavy jet mass p, the C-parameter, the wide and total 
jet broadenings B\y and Bt and the transition from three-jet to two-jet final states in 
the Durham jet algorithm I3. After reviewing the theoretical framework and fixing the 
notation in Section 2, we present results for the NNLO coefficients of the first five moments 
of the above-mentioned event shape distributions in Section 3. The numerical impact of 
these corrections and implications for the description of experimental event shape data are 
discussed in Section 4. 



2. Theoretical framework 

The perturbative expansion for the distribution of a generic event shape observable y up 
to NNLO at centre-of-mass energy yfs for renormalisation scale /i 2 = s and a s = a s (s) is 
given by 

1 da /a s \ dA /a,\ 2 dB /a s \ 3 d(7 



a had dy feJ d^ + feJ d^ + Kto) d^ + ° (aJ • (2 ' 1} 

Here the event shape distribution is normalised to the total hadronic cross section <7h a d- 
With the assumption of massless quarks, then at NNLO we have, 

^=-»( 1+ ^(£) + ^(£) 2+0 ^)) • < 2 ' 2 > 

where the Born cross section for e + e~ — > qq is 

a = —Ne 2 q . (2.3) 
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The constant K 2 is given by [47-49] 



3 2 

"jpF + CfCa 



/ _ _ 44C 3 j + CfTrNf (-22 + 16C 3 : 



where the QCD colour factors are 

C A = N, C F 



N 2 



2N 



Tr 



1 

2 



(2.4) 



(2.5) 



for N = 3 colours and Np light quark flavours. 

The perturbative QCD expansion of (y n ) is then given by 



V2vr 



2tt 



C y , n + 0(a s ) 4 



(2.6) 



In practice, we compute the perturbative coefficients An, B n and C n , which are all 
normalised to uq: 

"4y,n, £>y,n and C Vt n are straightforwardly related to A y ^ n , B y . n and by 

Ay^n = Ay^ n , 



By,n — By,n ^ ^ F "^V< n ' 
Cy,n = Cy,n ~ ^CpBy^ + [ -C F - K 2 ) Ay, n ■ 



(2.8) 



These coefficients are computed at a renormalisation scale fixed to the centre-of-mass en- 
ergy, and are therefore just dimensionless numbers for each observable and each value of 
n. 

The computation of the coefficients is carried out using the parton-level event gener- 
ator program EERAD3 [3,50], which contains the relevant matrix elements with up to five 
external partons [51-54], combined using an infrared antenna subtraction method [55]. A 
recently discovered inconsistency in the treatment of large-angle soft radiation terms [6] 
in the original EERAD3 implementation has been corrected. They account for an initially 
observed discrepancy between the EERAD3 results and the logarithmic contributions (com- 
puted within SCET) to the thrust distribution to NNLO [56], which are now in full agree- 
ment. 

In terms of the running coupling a s (fi 2 ), the NNLO expression for an event shape 
moment measured at centre-of-mass energy squared s becomes: 

2 



2vr 

a s (n) 



A 



+ < 2. 



y,n 

3 



2vr 



By,n + Ay, n Po log 



'y,n ' s-iy,; 



C y , n + 2B y , n Po log + Ay, n Po log —+Pl log 



.2 V 



(21 
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3. NNLO contributions 



The perturbative LO and NLO coefficients Ay >n , By,n have been known for many years [57- 
60], and have been used extensively in the experimental studies of LEP data. Our major 
new results presented here are the perturbative NNLO coefficients C y , n . These coefficients 
receive contributions from six different colour factors: 

N 2 , N° , l/N 2 , Np N , N F /N, N F . 

We have computed these six different contributions separately. 

The evaluation of the perturbative coefficients Ay ;n , B y>n and C y ^ n requires the numer- 
ical integration of the weighted perturbative event shape distributions down to the exact 
two-jet boundary y = 0. At this boundary, all event shape distributions diverge like 1/y, 
and are further enhanced by large logarithmic corrections at NLO and NNLO. Since the 
integrand for the n-th moment is the distribution weighted by y n , the moments themselves 
are finite. The first moment does however contain an integrable logarithmic singularity, 
and thus receives sizable contributions from the region y — > 0. 

The numerical convergence of the perturbative coefficients computed by EERAD3 dete- 
riorates for y — ► for the following reason. The four-particle and five-particle phase spaces 
in EERAD3 are generated with a lower cut-off So on all dimensionless invariants Sij/s. This 
cut-off is required to avoid sampling kinematical regions where the numerical evaluation 
of the four-parton one-loop matrix element becomes unstable due to the presence of large 
inverse Gram determinants, and also to avoid regions where the antenna subtraction pro- 
cedure induces cancellations over too many orders of magnitude for the remainder to be 
evaluated reliably. The independence of the results on So serves as a check for the proper 
implementation of the subtraction. In our calculation of the event shape distributions [1], 
we used So = 10~ 5 as a default, observing that this value allowed a reliable evaluation 
of all contributions while ensuring numerical stability and convergence. In approaching 
y — > So from above, the unresolved invariants are no longer small compared to the resolved 
invariants, and the calculation is likely to miss parts of the finite remainder from the sub- 
traction. For y < So, no reliable evaluation is possible. Consequently, we have to impose a 
lower cut yo on y when evaluating the perturbative coefficients. This cut is likely to affect 
especially the first moment. Due to the strong suppression of the y — > region in the 
higher moments, this cut does not play a role for n > 2. 

Therefore, in evaluating the perturbative coefficients, especially for the first moment, 
we carried out systematic studies to determine the values for the technical parameters So 
and yo- It is observed that yo = 10~ 5 is sufficient for all observables. In approaching 
this value, one observes already a reasonable convergence of all colour factor contributions 
to the first moments, and from studies with lower values of yo < 1CT 5 , we estimate the 
contribution from this region to each colour factor not to exceed five per cent of the total. 
The choice of So depends on the event shape variable under consideration. Already at 
NLO, one observes that So = yo/W is sufficiently small for a reliable evaluation of the 
first moment of (1 — T), p, C and I3, while By/ and Bt require much smaller So- In 
the distribution of both broadenings, one observes that a cut So induces a modification of 
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the distributions up to y ~ VSo, thereby rendering the numerical evaluation of the first 
moment of the broadenings much more challenging. This effect is illustrated in Figure 1 
for the NLO-contribution to the B\y distribution, evaluated for different values of <5o- For 
the NNLO contributions to the first moment of the broadenings, we chose 5q = 1CT 7 , 
which is clearly a compromise between the accuracy of the evaluation and the numerical 
convergence. We have verified for both broadenings that for fixed yo = 10 -5 , the numerical 
values of first moments do stabilize at 6q = 10~ 7 , and we estimate that the individual 
colour factor contributions do not change by more than three per cent for lower 5q. Full 
studies with much lower 5o were however not possible because of rapid deterioration of the 
numerical convergence. For these reasons, a systematic error of about 50% on the first 
moments of the broadenings should be taken into account. 

In Tables 1 and 2, we list the individual colour factor contributions to the first five 
moments of the different event shape variables. Table 3 contains the total LO, NLO and 
NNLO coefficients. The moments given in Tables 1 to 3 have been calculated with 5q = 1CP 7 
for Bw and Bt, where studies of other 5q values lead us to a systematic error estimate of 
10% (except for the first moment, where the systematic error is estimated at 50%). For 
the other event shape variables, the default for the first moments is Sq = 10 -6 with a 
systematic error of 6%, while the values for the higher moments have been calculated with 
5o = 10 -5 and studies of other 8o values result in a systematic error estimate of 3% for 
(1 - T), p, C and Y 3 . 

The behaviour in the different event shape variables is best seen by considering the 
ratios 

Knlo = ,4y ' w + ® V ' n fe) (3.1) 
and _ 

•Ay,n 

which we display for a s = 0.124 in Figures 2-4. Although substantially above the world 
average, it turned out that this value of a s is the best-fit result [42] of fitting the pure 
NNLO predictions to event shape data from ALEPH. Consequently, we also employ this 
value in the event shape studies here. We observe that the perturbative corrections affect 
the different shape variables in a substantially different manner, which was already evident 
at NLO (see for example [13]). While the size of the corrections is largely independent of n 
for p, Bw and F3 (and roughly K. ~ 1.3 at NNLO for all three variables), it increases with n 
for (1 — T), C and Bt (amounting to K, ~ 3 for n = 5 for all three variables). This picture is 
consistent with the observations made on the perturbative stability of the distributions [1], 
where the former three displayed uniform and moderate NNLO corrections, resulting in a 
stabilization of the perturbative expansion, while the latter three suffered large corrections. 
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4. Numerical results and comparison to data 

Moments of event shapes have been measured by various e + e~ collider experiments at 
centre-of-mass energies ranging from 10 GeV to 206 GeV [61]. These data are particularly 
interesting in view of the determination of non-perturbative power corrections. Since the 
power correction contributions scale like inverse powers of the centre-of-mass energy Q, 
one needs a large range in energies to empirically disentangle their contribution from the 
perturbative contribution, which scales logarithmically with the energy. 

However, even based on data covering about a factor twenty in energy range, missing 
higher order contributions may be misidentified as power corrections. The accuracy on 
the extraction of power correction contributions from data is therefore inherently linked to 
the precision of the perturbative prediction. We therefore expect the NNLO corrections 
computed here to lead to a considerable improvement of the power correction extraction 
from event shape moments. 

In Figure 5, we computed the energy dependence of the first moments of all event 
shape variables for as(Mz) = 0.124 at LO, NLO and NNLO. The theory uncertainty at 
each order is estimated by varying the renormalization scale in the strong coupling constant 
in the range Q/2 < [1r < 2Q. The theory predictions are compared to data from JADE 
and OPAL [61]. It can be clearly seen that the inclusion of the NNLO corrections leads 
to a substantial reduction (by a factor two to three compared to NLO) of the theoretical 
uncertainty in all variables. The theory error bands at NLO and NNLO generally overlap 
(although the overlap is only marginal for (1 — T), C and B?, which have the largest 
corrections, see above). We observe that the energy dependence of the first moment is 
well described by the purely perturbative NNLO QCD contribution for Y-j, while all other 
observables require the presence of additional non-perturbative contributions at low centre- 
of-mass energies. Comparing NLO and NNLO predictions, we observe that the NNLO 
contribution comes closer to the experimental data at low Q by about 30% for 1 — T and 
C, while only little improvement over NLO is observed for the other observables. The low- 
Q discrepancy between data and NNLO theory is most pronounced in p, C and (1 — T). It 
is especially interesting to note that the required size of the power correction contribution 
is uncorrelated with size of the perturbative higher order corrections: both p and I3 receive 
only moderate NLO and NNLO corrections; to describe the low-energy data on the first 
moment of p, one must postulate substantial (larger than in any other observable) power 
corrections, while data on I3 are well described without power corrections. 

The higher moments of the six standard event shapes are displayed in Figures 6 (1 — T), 
7 (p), 8 (C), 9 (Y 3 ), 11 (B T ) and 10 (B w ). The predictions are compared to JADE 
and OPAL data [61]. The qualitative behaviour of the higher moments of the different 
shape variables is similar to what was observed for the first moments. For (1 — T), C 
and Bt, we again observe that the NLO and NNLO theory error bands merely overlap. 
Despite the increase in the size of the NNLO corrections to the higher moments, this 
situation does not deteriorate for higher n, mainly because of the equally large increase 
of the NLO corrections. This behaviour should however be taken as an indication of 
the poor perturbative convergence of these three observables. For all n, we observe a 
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substantial discrepancy between NNLO theory and experimental data at low-Q in p, C and 
(1 — T) and Bt, indicating large power corrections to these observables. The relative size 
of the discrepancy does not decrease with n, illustrating the relevance of power corrections 
throughout the full kinematical range of the event shape variable (and not only in the 
two-jet region). Data on the higher moments of By/ and Y3 are well described by NNLO 
theory, leaving only little room for power correction contributions to these observables. 

With the new NNLO corrections, the theory uncertainty is reduced to the few per cent 
level for the moments of all event shape observables. It should therefore now be possible 
to perform precise studies of power corrections on these data sets. 

5. Conclusions and Outlook 

In this paper, we have computed the NNLO corrections to moments of the six most com- 
monly studied event shape variables: (1 — T), p, C, Bt, By/ and Y3. Especially the first 
moment receives substantial contributions from kinematical configurations deep inside the 
two-jet region, well beyond the region probed by event shape distributions. Consequently, 
the evaluation of the first moment is technically demanding, in particular for the jet broad- 
enings. 

The NNLO corrections to the first moment are moderate for all event shapes, leading 
to a /C-factor of order unity (in the range 0.95 to 1.4 for the different shapes). For higher 
moments, we observed that the NNLO corrections to the moments of (1 — T), C and Bt 
increase with increasing moment number n, resulting in large NNLO /C-factors of up to 
3. For By/, p and Y3, we observe that the corrections remain moderate towards higher n, 
indicating the better perturbative stability of these three variables. 

The energy dependence of the event shape moments can be used for the empirical 
determination of power corrections due to non-perturbative effects. For these studies, 
a reliable description of the perturbative contributions is mandatory. We observe that 
the newly computed NNLO corrections lead to a substantial reduction of the theoretical 
uncertainty on these predictions. From a comparison with data from JADE and OPAL, we 
observe that the energy dependence of the moments of By/ and Y3 can be largely described 
without any power corrections, while the other observables show clear indications for power- 
like contributions at low centre-of-mass energy. 

In view of the newly computed NNLO perturbative corrections, it will be very inter- 
esting to revisit the existing non-perturbative power correction models. Given the con- 
siderably reduced uncertainty on the perturbative prediction, a more precise extraction of 
hadronization parameters can now be envisaged. 
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<(i-m 


(P n ) 


(C n ) 






0?) 


N 2 n = 1 
2 
3 
4 
5 


2301 ± 11 
337.1 ±0.6 

63.2 ±0.2 
13.94 ± 0.02 
3.490 ± 0.005 


993 ± 12 
120.5 ±0.6 
15.7 ±0.1 
2.42 ± 0.03 
0.438 ± 0.006 


8816 ± 37 
4413 ± 5 
2381 ± 2 
1471 ± 1 
984 ±2 


348 ± 34 
174 ±1 
23.2 ±0.2 
3.20 ± 0.03 
0.456 ± 0.005 


3007 ± 59 

745 ± 1 
145.6 ±0.2 
32.13 ±0.03 
7.86 ±0.01 


470 ±5 
45.6 ±0.3 
6.64 ±0.07 
1.20 ±0.01 
0.249 ± 0.002 


iV° n = 1 
2 
3 
4 
5 


34.68 ±0.12 
-19.62 ±0.01 
-4.630 ± 0.002 
-1.1240 ±0.0006 
-0.2938 ±0.0001 


85.6 ±0.2 
-2.10 ±0.01 
-0.564 ± 0.003 
-0.0973 ± 0.0009 
-0.0163 ±0.0003 


285.4 ±0.4 
-220.1 ±0.2 
-165.13 ±0.06 
-113.92 ±0.03 
-80.72 ± 0.03 


4.0 ±4.4 
7.66 ± 0.02 
-0.028 ± 0.002 
-0.0474 ± 0.0007 
-0.2938 ± 0.0001 


383 ±2 
-21.64 ±0.03 
-8.656 ±0.005 
-2.285 ± 0.002 
-0.6025 ± 0.0005 


14.38 ± 0.06 
-1.181 ±0.005 
-0.243 ± 0.002 
-0.0494 ± 0.0005 
-0.0102 ±0.0001 


1/iV 2 n = l 
2 
3 
4 
5 


0.202 ± 0.008 
-0.2270 ± 0.0008 
-0.0010 ±0.0002 
0.00458 ± 0.00004 
0.00197 ±0.00001 


7.51 ±0.03 
-0.0799 ± 0.0009 
-0.0133 ±0.0003 
-0.00257 ± 0.00006 
-0.00065 ± 0.00002 


3.31 ±0.03 
-4.900 ± 0.007 
-0.931 ±0.006 
0.044 ± 0.004 
0.308 ± 0.003 


85.6 ±0.2 
0.232 ± 0.001 
-0.0089 ± 0.0002 
-0.00196 ± 0.00003 
-0.00036 ± 0.00001 


16.17 ±0.05 
-2.620 ± 0.004 
-0.3140 ± 0.0006 
-0.0467 ±0.0001 
-0.00790 ± 0.00003 


0.5788 ± 0.0041 
-0.0205 ± 0.0007 
-0.00432 ± 0.0001 
-0.00077 ± 0.00004 
-0.00019 ±0.00001 



Table 1: Individual colour factor contributions to event shape moments at NNLO: N 2 , TV , 1/N 2 . The given errors are the statistical errors from 
the numerical integration. Systematic errors originating from technical cut parameters (see text) are not shown in the table. 
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0T> 


N F N n = l 
2 
3 
4 
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-1716 ± 18 

-164 ± 1 
-26.1 ±0.2 
-5.23 ±0.04 
-1.202 ±0.008 


-1085 ± 15 
-95.8 ±0.9 
-12.8 ±0.2 
-2.30 ±0.04 
-0.47 ±0.01 


-6930 ± 80 
-2244 ± 11 
-1034 ±6 
-572 ± 4 
-351 ±3 


-664 ± 33 
-172 ±2 
-21.6 ±0.2 
-3.21 ±0.04 
-0.56 ±0.01 


-3279 ± 41 
-434 ± 2 

-70.6 ±0.3 
-14.06 ±0.06 
-3.23 ±0.02 


-451 ± 8 
-36.7 ±0.5 
-5.42 ±0.08 
-1.05 ±0.02 
-0.218 ±0.004 


N F /N n = l 
2 
3 
4 
5 


-6.67 ±0.01 
7.74 ±0.02 

1.478 ± 0.004 
0.3179 ± 0.0005 
0.0763 ± 0.0003 


-67.92 ±0.11 
1.39 ±0.02 

0.312 ± 0.003 
0.0596 ± 0.0006 
0.0129 ± 0.0002 


-85.48 ±0.07 
102.9 ±0.1 
59.29 ±0.09 
36.22 ± 0.05 
23.49 ± 0.03 


-465.2 ±0.3 
-3.16 ±0.02 
0.133 ±0.003 
0.0398 ± 0.0004 
0.0072 ± 0.0001 


-182.6 ±0.2 
23.79 ± 0.03 
4.903 ± 0.008 
1.075 ± 0.002 
0.2565 ± 0.0005 


-6.33 ±0.02 
0.642 ± 0.007 
0.123 ±0.001 
0.0241 ± 0.0004 
0.0052 ± 0.0002 


Np n=l 
2 
3 
4 
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253.9 ±0.2 
12.14 ±0.04 
1.388 ±0.007 
0.213 ±0.002 
0.0368 ± 0.0004 


271.32 ±0.19 
15.00 ± 0.04 
2.001 ± 0.008 
0.364 ± 0.001 
0.0760 ± 0.0004 


1123.5 ±0.5 
173.2 ±0.4 
55.8 ±0.2 
22.2 ±0.1 
9.15 ±0.07 


738.3 ± 0.3 
33.64 ±0.07 
3.829 ± 0.007 
0.609 ± 0.002 
0.1129 ±0.0004 


710.4 ±0.2 
27.79 ± 0.05 
2.42 ± 0.01 
0.240 ± 0.002 
0.0111 ±0.0005 


83.24 ± 0.05 
5.50 ±0.02 

0.792 ± 0.003 
0.1479 ±0.0008 
0.0320 ± 0.0003 



Table 2: Individual colour factor contributions to event shape moments at NNLO: Np N, Np/N, Np. The given errors are the statistical 
from the numerical integration. Systematic errors originating from technical cut parameters (see text) are not shown in the table. 
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A n = 1 
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2.1035 
0.1902 
0.02988 
0.005858 
0.001295 


2.1035 
0.1902 
0.02988 
0.005858 
0.001295 


8.6379 
2.4317 
1.0792 
0.5685 
0.3272 


4.0674 
0.3369 
0.04755 
0.008311 
0.001630 


4.0674 
0.3369 
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0.008311 
0.001630 


0.8942 
0.08141 
0.01285 
0.002523 
0.0005570 


£ n = 1 
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44.999 ± 0.002 
6.2595 ± 0.0004 
1.1284 ±0.0001 
0.24637 ± 0.00003 
0.06009 ± 0.00001 


23.342 ± 0.002 
3.0899 ± 0.0008 
0.4576 ± 0.0002 
0.08363 ± 0.00003 
0.01759 ± 0.00001 


172.778 ± 0.007 
81.184 ±0.005 
42.771 ± 0.003 
25.816 ± 0.002 
16.873 ± 0.001 


-9.888 ± 0.006 
4.5354 ± 0.0005 
0.6672 ± 0.0001 
0.10688 ±0.00002 
0.01865 ±0.00001 


63.976 ± 0.006 
14.719 ±0.001 
2.7646 ± 0.0003 
0.60690 ± 0.00006 
0.14713 ±0.00002 


12.689 ± 0.135 
1.2929 ±0.0001 
0.19901 ±0.00005 
0.03777 ± 0.00001 
0.00804 ± 0.00001 


C n = 1 
2 
3 
4 
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867.60 ±21.19 
172.80 ± 1.22 
35.38 ±0.26 
8.13 ± 0.04 
2.109 ±0.009 


204.55 ± 19.79 
38.93 ± 1.12 
4.626 ±0.196 
0.445 ± 0.051 

0.0367 ±0.0121 


3212.2 ±88.7 
2220.9 ± 12.0 
1296.6 ±6.7 
843.1 ±3.9 
585.0 ±3.2 


46.2 ± 47.8 
40.7 ± 1.9 
5.61 ±0.26 
0.591 ±0.051 
-0.273 ±0.011 


655.5 ± 72.0 
337.8 ±2.5 
73.34 ± 0.36 
17.06 ±0.07 
4.29 ±0.02 


110.6 ±9.3 
13.84 ±0.62 
1.892 ±0.107 
0.274 ± 0.027 
0.0585 ± 0.005 



Table 3: Contributions to event shape moments at LO, NLO, NNLO. 
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Figure 1: The NLO contribution to the wide jet broadening distribution Bw dB(Bw) / 'dBw eval- 
uated for different technical cut-offs 8q on all phase space variables. 
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Figure 2: The ratios Xnnlo (^nlo) of the NNLO (NLO) corrections relative to LO for the first 
five moments of the thrust and heavy mass distributions evaluated at /i = Q with a s = 0.124. 
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Figure 3: The ratios Xnnlo (-^nlo) of the NNLO (NLO) corrections relative to LO for the first 
five moments of the C-parameter and total jet broadening distributions evaluated at /i = Q with 
a s = 0.124. 
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Figure 4: The ratios Xnnlo (-^nlo) of the NNLO (NLO) corrections relative to LO for the first 
five moments of the wide jet broadening and Y3 distributions evaluated at fi — Q with a s = 0.124. 
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Figure 5: First moments of six event shape variables. The bands show the scale variations obtained 
by using // = 2Q and /i = Q/2. The data are from the JADE and OPAL experiments, taken 
from [61]. 
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Figure 6: Perturbative corrections to the higher moments of 1 — T . The bands show the scale 
variations obtained by using fi = 2Q and /i = Q/2. The data are taken from [61]. 
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Figure 7: Perturbative corrections to the higher moments of the normalised heavy jet mass p. 
The bands show the scale variations obtained by using [i — 2Q and p — Q/2. The data are taken 
from [61], which only contains data for the moments of Mh = ^fp-, such that no data are shown for 
n = 3,4,5. 
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Figure 8: Perturbative corrections to the higher moments of the C-parameter. The bands show 
the scale variations obtained by using /i = 2Q and fj, = Q/2. The data are taken from [61]. 
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Figure 9: Perturbative corrections to the higher moments of the 2-jet to 3-jet transition parameter 
I3. The bands show the scale variations obtained by using fj, — 2Q and /i = Q/2. The data are 
taken from [61]. 
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Figure 10: Perturbative corrections to the higher moments of the wide jet broadening By/ . The 
bands show the scale variations obtained by using fi = 2Q and /i = Q/2. The data are taken 
from [61]. 
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Figure 11: Perturbative corrections to the higher moments of the total jet broadening Bt- The 
bands show the scale variations obtained by using /i = 2Q and fi = Q/2. The data are taken 
from [61]. 
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